[Analog_Engine.git] / scripts / TP1 08.22 Hollow Cylinder.AESL

IDENTIFICATION DIVISION
PROGRAM-ID HollowCylinder
VERSION 20240312
COMMENT A homogeneous hollow cylinder (mass m, ration inner to outer radius epsilon) is connected to a massless spring (constant lambda). 
COMMENT The system is deflected from its rest position by the distance a. 
COMMENT The cylinder rolls without sliding. 
COMMENT x'' = - 2*lambda * 1/m * 1/(3+epsilon²) * x

ENVIRONMENT DIVISION
ENGINE Anabrid-THAT
TIMEBASE 1ms
REQUIRES COEFFICIENT 4, INTEGRATOR 2, INVERTER 1

DATA DIVISION
OUTPUT OUTPUT.X x
COEFFICIENT.1 2*lambda
COEFFICIENT.2 1/m
COEFFICIENT.3 1/(3+epsilon²) # = 1/3 (0,333) for epsilon=0 and 1/4 (0,25) for epsiolon=1 or anything in between
COEFFICIENT.4 a 

PROGRAM DIVISION
+1 -> COEFFICIENT.A -> a

x'' -> INTEGRATOR -> -x'
-x', IC:a -> INTEGRATOR -> x
x -> COEFFICIENT.2*lambda -> 2*lambda*x
2*lambda*x -> COEFFICIENT.1/m -> 2*lambda*1/m*x
2*lambda*1/m*x -> COEFFICIENT.1/(3+epsilon²) -> 2*lambda*1/m*1/(3+epsilon²)*x
2*lambda*1/m*1/(3+epsilon²)*x -> INVERTER -> -2*lambda*1/m*1/(3+epsilon²)*x = x''

OPERATION DIVISION
MODE REPEAT
OP-TIME 52ms
Commit	Line	Data
6e4089e8 P	1	IDENTIFICATION DIVISION
	2	PROGRAM-ID HollowCylinder
	3	VERSION 20240312
	4	COMMENT A homogeneous hollow cylinder (mass m, ration inner to outer radius epsilon) is connected to a massless spring (constant lambda).
	5	COMMENT The system is deflected from its rest position by the distance a.
	6	COMMENT The cylinder rolls without sliding.
	7	COMMENT x'' = - 2lambda 1/m * 1/(3+epsilon²) * x
	8
	9	ENVIRONMENT DIVISION
	10	ENGINE Anabrid-THAT
	11	TIMEBASE 1ms
	12	REQUIRES COEFFICIENT 4, INTEGRATOR 2, INVERTER 1
	13
	14	DATA DIVISION
	15	OUTPUT OUTPUT.X x
	16	COEFFICIENT.1 2*lambda
	17	COEFFICIENT.2 1/m
	18	COEFFICIENT.3 1/(3+epsilon²) # = 1/3 (0,333) for epsilon=0 and 1/4 (0,25) for epsiolon=1 or anything in between
	19	COEFFICIENT.4 a
	20
	21	PROGRAM DIVISION
	22	+1 -> COEFFICIENT.A -> a
	23
	24	x'' -> INTEGRATOR -> -x'
	25	-x', IC:a -> INTEGRATOR -> x
	26	x -> COEFFICIENT.2lambda -> 2lambda*x
	27	2lambdax -> COEFFICIENT.1/m -> 2lambda1/m*x
	28	2lambda1/mx -> COEFFICIENT.1/(3+epsilon²) -> 2lambda1/m1/(3+epsilon²)*x
	29	2lambda1/m1/(3+epsilon²)x -> INVERTER -> -2lambda1/m1/(3+epsilon²)x = x''
	30
	31	OPERATION DIVISION
	32	MODE REPEAT
	33	OP-TIME 52ms